Electron space charge effect on spin injection into semiconductors 
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We consider spin polarized transport in a ferromagnet-insulator/semiconductor/insulator- 
ferromagnet (F1-I-S-I-F2) junction. We find that the spin current is strongly dependent on the 
spin configurations, the doping and space charge distribution in the semiconductor. When the 
ferromagnet-semiconductor interface resistance is comparable to the semiconductor resistance, the 
magnetoresistance ratio of this junction can be greatly enhanced under appropriate doping when 
the space charge effect in the nonequilibrium transport processes is taken into consideration. 



PACS numbers:73.40.-c,71.70.Ej,75.25.+z 

There is much interest recently in incorporating mag- 
netic elements into semiconductor structures. One of the 
current focus lies in injecting spin polarized electrons into 
non-magnetic semiconductors [1—7]. This is partially mo- 
tivated by the high magnetoresistance observed in ferro- 
magnet tunnel junctions [8]. However, experiments have 
so far observed a small magnetoresistance ratio of 1% 
[9,10] in F-S-F structures. Schmidt et al [11] emphasized 
that the large resistance of the nonmagnetic semiconduc- 
tor dominates the total resistance and thus lowers the 
magnetoresistance ratio. A similar discussion based on 
a ballistic picture was given by Grundler [12] and Hu 
et al [7]. Rashba pointed out that in a ferromagnetic 
metal/insulator/semiconductor (FIS) junction the resis- 
tance of the barrier insulator layer can be comparable 
to that of the semiconductor, AR/R may be raised by a 
considerable percentage [13]. Many different geometries 
of the tunneling junction were discussed in Ref. [14]. 

Although it is possible that to lift the magnetoresis- 
tance to near 10% with a very high electric field [15], 
there is no a low field scenario yet. On the other hand, the 
effect of the Coulomb interaction between the electrons 
in the non-magnetic semiconductor in the nonequilibrium 
transport process was not taken into consideration in the 
above mentioned works. In previous works on the ferro- 
magnetic metal/insulator /ferromagnetic metal junction 
[16], one of us found that this interaction between the 
electrons in nonequilibrium transport processes played an 
important role in the characteristics of the device. Un- 
der steady state nonequilibrium conditions, a magnetiza- 
tion dipole layer much larger than the charge dipole layer 
is induced at the interface. The magnetization dipole 
layer is zero under equilibrium conditions. In this paper, 
we consider the role played by the Coulomb interaction 
in an F1-I1-S-I2-F2 junction. To completely solve the 
nonequilibrium transport processes with interactions is 
highly non-trivial. It requires solving self consistently 
Boltzmann type spin transport equation with the Pois- 
son equation. What we want to touch in this work is 
using a simple Hartree approximation to check if the ef- 
fect comes from the interactions is significant. We found 
that the spin polarized current is strongly affected by the 
interaction in this approximation. Because of the large 



difference in the spin-dependent current between the par- 
allel and anti-parallel configurations, which is caused by 
the space charge effect, the total resistance becomes de- 
pendent on the spin configuration. We found that mag- 
netoresistance ratio AR/R can be increased by one order 
of magnitude higher than the result obtained without the 
interactions under appropriate doping of the semiconduc- 
tor. Typically, AR/R ~ 3% for the non-interaction elec- 
trons while AR/R may attain values as high as ~ 50% 
due to the space charge. (See, Figures 2 and 3 below.) 
We also study another physical observable, the polariza- 
tion P = 3 j + + j~ an d find that it is strongly dependent 
on the space charge distribution. We now describe our 
result in detail. 

General description: We consider the junction Fl-Il-S- 
I2-F2. The thickness of the metal (Fl, F2), the insulating 
barriers (II, 12) and the semiconductor (S) are denoted 
by L L ' R , di.2 and x , respectively (See, Figure 1). For 
the practical case, Xq is less than the spin diffusion length 
in the semiconductor In- The charge screening lengths 
in the metal and the semiconductor are denoted by Xl.r 
and Xn, respectively. In our model, we assume A at <C 
In and \r,l Ir,Li the spin diffusion lengths in the 



metal. Typically, A 



io- 



l ~ 10 nm, and In > 



1/^m. xq ~ lOOnrn to lpm, depending on the structure 
of junctions. The screening length in the semiconductor, 
A at, is dependent on the doping of the semiconductor and 
can vary in a wide range, say lOnm for the heavy doped 
semiconductor and 100nm-l/xm for the lightly doped or 
undoped semiconductor. 

The problem we would like to solve is defined by the 
following four sets of equations. 

(1) The first is the total charge-current conservation, 



v-j = --^ 



dp 

at' 



(i) 



where p is the charge density. 

(2) The second is the diffusion equations of the spin s- 
dependent current. In order to allow an analytic analysis, 
we take a simple Hartree approximation into account to 
see the space charge effect. Under such an approxima- 
tion, the diffusion equations reads [16] 
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j s = a s (Vp s -VW + E), 



(2) 



where the magnitude of the electric charge has been 
set as one; E is the external electric field; the chemi- 
cal potential p s is related to the charge density p s by 
V/k s = -^p- where N s is the spin-dependent density of 
state. Wo = / dr , Ui n t{r — r')p{f l ) is the potential caused 
by the screened Coulomb interaction Ui nt (r) and V Wo is 
the so-called screening field induced by Ui n t(f). (In the 
present device geometry, the three dimensional integra- 
tion is reduced to a one-dimensional one.) 
(3) The third is the magnetization relaxation equation 
where the magnetization density M = pj — pi relaxes 
with a renormalized spin diffusion length I. In the relax- 
ation time approximation, one has 



W 2 M-M/l 2 = 0. 



(3) 



(4) The fourth is the boundary condition at the interfaces 

A/2« - AW = r(l - sj)j s , (4) 

where p s = p s + Ex; Ex is the voltage drop on the left 
side of the barrier and AW is the electric potential drop 
across the barrier, which is assumed much smaller than 
A/i s ; r s = r(l — sj) is the barrier resistance. We assume 
that there is no spin relaxation in the insulator, the spin- 
dependent currents are continuous across the junctions, 
jH-di/2) = j^(d 1 /2) and jf (x + d 2 /2) = j?(x - 
d 2 /2). 

In addition, we have the neutrality condition for the 
total charges (Ql,n,r, e.g., Qn — Jd"/^ 1 ^ pdx) accu- 
mulated at the interfaces. By Gauss' law, for the point 
d\/2 < x < x + di/2, i.e., inside the semiconductor, the 
potential Wq is determined by 



VWofi) = 4ttQ l + 4tt 



pdx, 



di/2 



(5) 



whose constant part of the right hand side gives the con- 
straint on the charge while the rc-dependent part gives the 
function form of the potential Wo. Another constraint is 
the neutrality of the system: 



Ql + Qn + Qr = 0. 



(6) 



With those sets of equations (Eqs.(l)-(4)) and the two 
constraint on the charges (Eqs.(5) and (6)), the problem 
can be solved. The formal solutions of the problem are 



p L (x) = -k p L oe (x+dl/2)/\L ::± n L (,■ ,/, 2) /, 



A 2 



(7) 



with similar solutions for the right hand side. In the 
semiconductor, if Ajy «iOi 



p W {x) = ^ p W e -(-^/2)/A» + ^L^) e -(x-«J 1 /2)/I W) 

In l N 

pM {x) = *N p (2) e ( X - X0+d2 /2)/X> 

In 



, V (2) (x-x +d 2 /2)/l N 
iV 



(8) 



can be obtained similarly. All of coefficients in 
(7) and (8) can be determined by using Eqs.(l)-(4) and 
the constraint (5) and (6). The screening potential Wo 
is determined by Gauss' law. The total current is j = 
^2 s j s - Although j s are not a constant, the total current 
j is still a constant. 

The spin- dependent current: To demonstrate the essen- 
tial physics, we simplify matters by setting the parame- 
ters of the metals and barrier widths on the left and right 
sides to be the same: Xr, = Xl = A, II = Ir = I, d\ = 
d 2 = d and so on. The resistances of the barrier layers 
arc taken as = r^ = r; 71 = 72 = 7 for the parallel 
configuration and 71 = —72 = 7 for the anti-parallel con- 
figuration. To illustrate, we focus on the calculation in 
the left barrier located at x — 0. The tunneling resistance 
is given by r ( p = r(l - 7s) = r , s exp[d(/c a (//) - k s (0))], 
where K s (p) oc J Q dx[2m(U — A / u s (0)a;/d)] 1//2 , with U 
the barrier height. That is K s {p) = no,s gA^Tol ft — — 
A/i s (0)) 3 / 2 ]. Here Ap s = Ap s /U. The current j^(x) at 
x = is dependent on the bias voltage and the space 
charge distribution, which is given by 



where 



1 + 



j s L (0) = A s j 0s 



4ttX 2 a{0 - s) pfb + Mq 



I l-a[3 



E 



(9) 



(10) 



and jos = a s E is the current with no interaction; (3 (a) 
measures the spin asymmetry of the conductivities a s 
(densities of states at the Fermi surface N s ): <r s = ^(1 + 
ps) and N s = ^N F (1 + sa) where N F is related to the 
screening length A by = 2nX 2 \Z^p ■ Noting that both 
Pi and Mq are proportional to the external electric field 
E, A s is solely determined by the material parameters. 

Eq. (9) implies that the spin-dependent current j s (0) 
passing the interface differs a factor A s from the non- 
interacting current jo s (0) which, for the parallel configu- 
ration, is given by 



Jos(0) 



V 



RN + 2r .sY s (0)' 



(11) 



where Y s (0) = ^.^(i-Ci-A^o))^)-!]. For the 
anti-parallel configuration, 



V 



(12) 



2 



which is almost spin- independent. 

The magnetoresistance: Since the total current is con- 
stant everywhere, we have, for the parallel configuration 

i-EJ.w/v-EsTT&oj- (13) 

and for the anti-parallel configuration, we have 

^^'■'■'^ ^rsb-i ' (14) 

From these, we obtain the magnetoresistance ratio 
AR = R AP - R P 2X 
R ~ (Rap + Rp)/2 2 + X : [ ' 

^ Af(R N + Zs'ro,s'Y s t p (0)) 
4- 2(R N + 2r , s Yf(0)) 

For non-interacting electrons, = A p = 1 and one can 
see that AR/R will not be beyond a maximal value (at 
V -> 0) about 3.23% for r 0+ : r _ : R N = 1 : 2 : 1 and 
decays as the bias voltage increases. In Fig. 2, we plot 
the bias-dependence of the magnetoresistance ratio for 
both the non-interacting and the interacting case with 
x ~ 1.25/xm. Here we take A = O.lnm, I = 20nm, 
a = f3 = 1/2 and In = 3/im [17]. After the space 
charge effect is included, the ratio (15) grows as Xn be- 
comes smaller. This screening length can be controlled 
by doping the semiconductor (or by applying an external 
gate voltage). The small An requires heavy doping of 
the semiconductor. For a heavily doped semiconductor, 
it may arrive at the giant magnetoresistance value, say 
56% as shown in the figure, which raises orders of the 
magnitude comparing to the magnetoresistance for the 
non-interacting electrons. 

There are some points to be emphasized in our calcu- 
lations, i) In the large Rn limit, AR/R tends to zero 
no matter if there is any interaction effect, ii) At zero 
bias, the magnetoresistance is highest. In Fig. 2, we 
show the dependence of the magnetoresistance ratio on 
Ajv/^o for different xq(^$> I) at zero bias [18]. For a given 
xo, we see the magnetoresistance is also dependent on I. 
iii) The magnetoresistance is reduced as the bias voltage 
increases. We expect that the effective barrier height U, 
say several voltage, for the FIS interface to be larger than 
that for the F-I-F structure because the work function of 
the semiconductor is smaller than that of a metal. So on 
an absolute scale of voltage, this effect is not as big as 
in the conventional F-I-F structure. Namely, the large 
magnetoresistance can be anticipated at an experimental 
relevant bias voltage. 

The polarization: We consider another experimental 
measurable quantity [5], the polarization of current, 
which is defined by P — J - + ^T'wn! ■ We first calcu- 

■' J+(0)-h7-(0) 

late P by ignoring the space charge effect. In the anti- 
parallel configuration P is nearly zero. The polariza- 
tion P for the parallel configuration is also small if Rn 



dominates but it may arrive a considerable magnitude 
if the ros is of the same order of Rn, e.g., if the ratio 
r o+ ■ ro- ■ Rn = 1 : 2 : 1, it may arrive at 40% at the 
zero bias. (In a wide range, our result is not very sensi- 
tive to this ratio if these resistances have the same order 
of magnitude.) P decays as the bias voltage increases, 
say for V/U <~ 1, P drops to less than 10%. 

To see the space charge effect, we look at A s . We found 
that A s are nearly independent of In if In is the longest 
length scale. We fix In — 3^m and n 0s d = 300. For 
the anti-parallel configuration, the polarization is given 
by P ap = , which is bias-independent. In a large 

screening length Xn and x , c.g, lOOnm and 1.5 yum, P ap 
is of the same order of the magnitude as that in the non- 
interacting case. However, for a small Xn and xq, e.g., 30 
nm and 0.5/zm, P ap can be larger than one because j + , 
j- can go in opposite directions. P ap becomes significant 
when Xn and xo are getting smaller. The physics behind 
this phenomenon can not be well-understood yet at this 
stage. An intuitive consideration is as follows: The ef- 
fective driving potential A/i s is different in sign between 
spin up and spin down electrons, because the net charge 
density at the interface is less than the magnetization 
density [16]. This inequality can be traced to stem from 
the fact that the magnetic length is much larger than the 
screening length. Because of this sign difference of the 
driving potential, the currents for the two spin compo- 
nents may go in opposite directions. However, this con- 
sideration has to be examined thoroughly and we shall 
leave it to the further works. The polarization of the par- 
allel configuration is bias dependent. For a large Ajv and 
small xo (say 100 nm and 0.5/xm), P p decays as the bias 
voltage increases. However, for a small Xn (say 30 nm) 
and large xo (say 2/zm), P p increases as the bias voltage 
and finally saturates. However, these bias dependent ef- 
fects happen only in the high bias regime, say V/U > 1. 
Thus, for the practical usage, this bias dependence may 
not be so important. 

In conclusions, we have considered the space charge 
effect in a multiple FIS junction. In the Hartrce approx- 
imation, we found a conspicuous difference of the spin- 
dependent current in the parallel configuration from that 
in the non-interaction models while the difference is mi- 
nor for the anti-parallel configuration. This leads to an 
order of magnitude growing in the magnetoresistance ra- 
tio in such an approximation from the non-interaction 
model if the parameters of the junction are appropri- 
ately chosen. Although our theory is based on a spe- 
cific junction geometry and a simple approximation, it 
has indicated that the electron interaction may play an 
important role in noncquilibrium transport processes for 
some types of spintronic devices on the spin injection. To 
quantitatively solve the nonoequilibrium transport pro- 
cesses, a stricter treatment is required. 
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Fig. 2 The magnetoresistances versus the bias voltage 
for non-interacting junction and the junctions with the 
space charge distribution. K 0s d = 300 has been taken. 

Fig. 3 The magnetoresistance versus Xn/xq in V = for 
given semiconductor sizes x . The curves corresponding 
to different I are plotted in each diagram. 
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Fig. 1 The sketch of the geometry of the junction. The 
quantities of the left (right) ferromagnet is indexed as 'L' 
('R') . 
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